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This  Materials  Engineering  for  Affordable  New  Systems  (MEANS)  program  focused  on 
developing  an  integrated  computational  framework  for  predicting  the  conditions  under  which 
turbine  blades  can  be  produced  that  are  free  of  deleterious  solidification  defects.  To  achieve 
this  ambitious  goal  we  studied  the  solidification  process  over  length  scales  from 
subnanometer  to  millimeter  and  produced  a  fundamental  experimental  and  theoretical 
framework  for  defect  prediction.  The  present  work  aimed  at  producing  a  fundamental 
predictive  theory  in  the  presence  of  these  effects  and  testing  this  criterion  experimentally.  A 
criterion  for  the  onset  of  channels,  freckles,  and  freckle  chains  in  nickel-base  superalloys 
relevant  to  the  aerospace  industry  has  been  developed  by  Davis  and  Voorhees.  This  criterion 
involves  the  Rayleigh  number  that  measures  buoyancy  in  the  mushy  layer  versus  the 
dissipative  effects  of  heat  conduction  and  viscous  stresses,  as  well  as  a  second  parameter  the 
Biot  number  that  measures  the  lateral  heat  loss  and  hence  the  curvature  of  the  interface.  The 
segregation  coefficients  and  density  parameters  of  multicomponent  supcralloys  were 
calculated  using  hybrid  a/>-/w7/o/eomputational-thermodynamics  modeling  by  Asta  and 
Woodward.  The  permeability  of  dendritic  arrays  in  superalloys  has  been  determined  using 
three-dimensional  reconstructions  of  the  solid-liquid  mush  and  finite-element  fluid 
simulations  by  Pollock  and  Spowart. 

Close  interaction  with  industry  ensured  that  computational  models  were  developed  in  a  form 
that  is  usable  in  their  efforts  to  design  new  alloys  and  processing  routes.  Given  the  broad 
research  focus  of  the  project  the  MEANS  team  was  composed  of  engineers  and  scientists 
from  government  laboratories,  industry,  and  universities  and  is  diverse  in  its  research 
expertise. 
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Background 


Preventing  the  formation  of  freckles  and/or  misoriented  grains  is  important  to  the 
manufacturing  and  performance  of  single  erystal  superalloy  components.  These  defeets  are 
the  result  of  solute-indueed  eonveetive  instabilities  that  oeeur  at  the  solid-liquid  interface 
during  direetional  solidification.  The  two  most  eommon  defeets  observed  are  isolated, 
individual  high-angle  mis-oriented  grains  and  freekle-ehains.  Although  a  number  of  faetors 
have  been  identified  as  contributors  to  the  formation  of  these  defeets,  ineluding;  eooling 
rates,  refractory  alloy  eontent  and  easting  size  and  geometry,  eonveetive  flow  during 
solidification  is  regarded  as  the  precursor  event.  As  a  result,  the  most  widely  used  criteria  for 
the  prediction  of  these  defeets  generally  consider  the  ratio  of  the  buoyant  to  frietional  forees, 
an  interaction  ehiefly  quantified  by  the  Rayleigh  number.  Such  a  prediction  requires  not  only 
knowledge  of  the  fluid-flow  conditions  in  the  melt  but  also  a  detailed  understanding  of  the 
geometrical  domain  in  whieh  the  melt  flows.  Unfortunately,  a  definitive  understanding  of 
fluid  flow  at  the  seale  of  direetionally  solidified  dendritie  structures  is  still  laeking.  This  is 
partly  due  to  the  many  faetors  governing  flow  at  this  level  that  arc  both  dynamie  and  difficult 
to  quantify.  These  faetors,  often  eaptured  in  the  Rayleigh  parameter  Ra,  Eq.  (1),  include 
composition,  segregation  behavior,  permeability,  and  flow  ehannel  geometry: 


av 


where  Ap/po  is  the  density  gradient  in  the  liquid,  g  is  acceleration  due  to  gravity,  K  is 
permeability,  L  is  the  height  of  the  mushy  zone  or  relative  length  seale  and  av  is  the  product 
of  thermal  diffusivity  and  kinematic  viseosity.  Alloy  chemistry  influences  the  Rayleigh 
number  by  means  of  its  inherent  effeet  on  the  density  gradient,  diffusivity  and  viseosity. 

This  Rayleigh  number  criterion  is  derived  from  a  theory  of  two-eomponent  solidification  and 
is  successfully  used  in  multi-eomponent  superalloys  if  the  thermophysieal  constants  for 
segregation  (and  henee  buoyaney)  and  the  permeability  are  known.  These  predictions  are 
good  when  the  mush-melt  front  is  horizontal;  when  it  is  not,  as  is  often  the  ease  in  praetiee, 
the  theory  is  non-predietive.  When  real  eastings  are  examined,  it  is  noted  that  ehannel  and 
henee  freekles  and  freekle  ehains  preferentially  form  near  the  mold  surfaee  where  heat  losses, 
interface  eurvaturc,  and  horizontal  permeability  variations  arc  greatest.  The  present  work  was 
aimed  at  producing  a  fundamental  predietive  theory  in  the  presence  of  these  effects  and 
testing  this  criterion  experimentally 

Densities  of  Liquid  Superalloys 

In  direetional  solidification,  the  mass-density  difference  between  the  hot  liquid  near  the 
solid-liquid  interface  and  the  eooler  melt  at  the  top  and  side  walls  of  the  easting,  is  the 
driving  foree  for  eonveetion  and  associated  formation  of  solidification  defeets.  The  mass- 
density  difference  thus  originates  from  both  the  composition  and  temperature  dependencies 
of  the  liquid-phase  molar  volume  (V(e,T)).  We  have  employed  ab-initio  moleeular  dynamics 
simulations  extensively  as  a  computational  framework  for  aiding  the  development  and 
validation  of  aeeurate  models  for  V(e,T)  in  Ni-based  superalloys. 
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AIMD  simulations  have  been  performed  to  derive  liquid-phase  atomic  volumes  for  twelve 
elemental,  binary  and  ternary  alloy  compositions  of  molten  Ni  with  Al,  W,  Ta,  and  Re 
additions  to  compute  atomic  volumes  at  temperatures  of  1830  and  1750  K.  Comparisons 
between  AIMD  results  and  direct  measurements  for  Ni,  Ni-Al,  Ni-Ta,  and  Ni-W  yield 
agreement  at  the  level  of  0.6-1 .8%,  with  the  calculated  volumes  being  systematically  larger 
than  experimental  volumes.  Simulations  were  performed  for  ternary  alloy  compositions 
containing  W,  Rc  and  Ta,  from  which  accuracy  limits  were  assessed  for  the  most  recently 
published  density  parameterization  at  alloy  compositions  outside  the  range  where 
experimental  data  is  presently  available. 

Overall  our  results  demonstrate  the  high 
accuracy  attainable  with  AIMD  and  the 
utility  of  this  method  for  providing  data 
where  experimental  results  are  unavailable. 

They  form  the  basis  for  refining  the 
development  of  density  models  for 
supcralloy  solidification.  In  addition,  the 
AIMD  results  present  a  wealth  of 
information  related  to  the  structure  and 
dynamics  of  molten  alloys,  useful  in  other 
contexts  related  to  the  processing  of 
supcralloys.  For  example,  the  calculated 
tracer  diffusion  coefficient  for  W  in  Ni-W 
molten  alloys  was  calculated  by  AIMD  to  be 
2.7(2)xl0'?  cm2/s,  in  excellent  agreement 
with  a  recently  published  value  of 
2.4(2)xl0‘5  cm2/s.  Further,  an  analysis  of 
the  structures  of  the  simulated  binary  Ni- 
based  alloys  shows  a  pronounced  tendency 
towards  chemical  short-range  ordering,  and 
a  tendency  towards  the  preferential 
formation  of  icosahedral  toplogical  short- 
range  order. 

In  collaboration  with  Jim  Lill  and  Chris  Woodward  at  AFRL,  efforts  in  the  latter  part  of  this 
project  led  to  the  implementation  constant-pressure  molecular  dynamics  into  the  commercial 
VASP  AIMD  code.  The  new  algorithms  implemented  in  this  work  laid  the  foundation  for 
extending  the  simulation  methodology  to  real  multicomponent  superalloy  compositions.  As 
an  example,  the  figure  above  shows  a  snapshot  from  a  simulation  of  an  eight-component 
model  RENE-N4  superalloy  which  formed  the  focus  of  the  experimental  effort  aimed  at 
validation  of  Rayleigh-number  criterion  for  the  prediction  of  the  formation  of  freckles  and 
related  solidification  defects. 


Figure  1:  Schematic  of  the  500  atom  cell 
(left)  and  the  evolution  of  the  internal 
pressure  (right)  over  a  2.6  psec  interval  for 
an  NPT  simulation  of  an  eight  component 
RENE-N4  alloy.  The  colors  of  the  atoms 
correspond  to  the  following  elements  Ni- 
light  grey,  Cr-green,  Co-gold,  Mo-purple, 
W-blue,  Al-red,  Ti-cyan,  Ta-black. 
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Permeability  via  Three-Dimensional  Reconstructions 

The  permeability  of  a  mush  is  a  strong  function  of  the  morphology  of  the  dendrites  and  the 
volume  fraction  of  liquid.  Because  3-D  imaging  of  dendrites  following  complete 
solidification  is  a  major  challenge,  our  approach  has  been  to  utilize  crystals  where 
solidification  has  been  abruptly  terminated  and  the  liquid  at  the  solidification  front  is  quickly 
decanted.  Preliminary  reconstructions  have  utilized  a  decanted  sample  of  single  crystal  Rene 
N4  supcralloy  solidified  in  the  University  of  Michigan  Bridgman  furnace.  Using  the 
prototype  RoboMET3D  system  developed  by  J.  Spowart  and  H.  Mullens  at  Wright-Patterson 
Air  Force  Base  in  Dayton  OH,  a  serial  sectioning 
technique  was  developed  for  the  Rene  N4  alloy. 

Optimum  slice  rates  of  approximately 
7sl  iccs/hour  were  reached  and  datasets  obtained 
were  a  collection  of  8-montagcd  images  taken  at 
10X  magnification  each.  The  resolution  of  each 
dataset  is  ,52um/pixel  and  in  total  the  size  of  each 
dataset  “slice”  was  on  the  order  of  153MB.  The 
recession  rate  between  each  imaging  set  was  on 
average  2.2um.  See  Figure  2  for  a  reconstructed 
image. 

The  three-dimensional  arrays  were  cropped  for 
regions  of  liquid  connectivity  and  surfaee  meshes 
were  generated  using  MIMICS  software  by 
Materialise  followed  by  volume  meshes  that  were 
created  in  GAMBIT  by  ANSYS.  After 
satisfactory  meshing  of  the  interdendritic  liquid 
was  accomplished,  global  smoothing  of  the  structures  was  introduced  to  dimmish  the 
pronounced  tertiary  dendritic  features  and  limit  flow  interaction  to  the  primary  and  secondary 
dendrite  arms  to  ensure  convergence.  FLUENT  was  used  to  simulate  fluid  flow  through  the 
dendritie  structure  in  the  following  manner.  A  direction  of  flow  was  assumed  and  a  pressure 
gradient  is  imposed  on  the  structure.  Boundary  conditions  of  zero  pressure  at  the  outlet  and  a 
flow  velocity  of  lOOpm/see  at  the  inlet  with  ‘no-slip’  at  the  walls  was  assumed,  and 
interdendritic  fluid  flow  was  assumed  to  be  steady  state.  Using  Darey’s  Law  the  permeability 
in  these  structures  were  then  calculated. 

For  eonvcctive  instabilities  to  develop  in  a  dendritic  structure,  flow  across  the  dendritie  array 
normal  to  the  solidification  direction  must  feed  the  plumes  that  flow  parallel  to  the 
solidification  direction  and  ultimately  result  in  freckles.  Thus  the  influence  of  dendritic 
structure  on  this  “cross-flow”  normal  to  the  solidification  direction  is  of  interest.  A  eross- 
sectional  volume  nominally  200  x  1000  x  1500  mm  was  selected  due  to  its  location  at  the 
dendrite  tips  as  well  its  high  connectivity  of  interdendritic  regions.  Streamlines  illustrating 
the  flow  directions  in  this  volumetrie  eross-seetion  are  shown,  Figure  3.  In  these  simulations, 
flow  is  assumed  as  steady  state  and  color  bar  illustrates  flow  speeds  throughout  the  mushy 
zone.  It  can  be  observed  that  flow  speed  increases  in  narrow  channels  where  flow  is 
constricted  and  the  flow  rate  must  increase.  Flow  also  seems  to  be  dominated  by  the  primary 


Figure  2:  Reconstructed  dendritic  solid- 
liquid  interface  of  Rene  N4,  gray 


represents  solidified  superalloy,  the  voids 
represent  interdendritic  regions  occupied 
by  liquid  prior  to  decanting. 


4 


and  secondary  dendrite  arms  found  in  the 
central  portions  of  the  region.  The  pertinent 
range  of  speed  magnitudes  observed  range  from 
175  pm/s  to  3500  pm/s.  These  values  appear 
reasonable  given  the  imposed  inlet  speed  of  100 
pm/s.  We  have  determined  the  permeability  of 
this  structure  to  be  1 . 1 6  x  1  O'  m  . 


Through  this  work  we  showed  that  the  primary 
and  secondary  dendritic  arm  spacings  of  the 
reconstruction  arc  consistent  with  two- 
dimensional  measures  and  are  within  the  range 
of  expected  values  based  upon  withdrawal  rates 
and  thermal  gradients.  In  addition,  for  the 
commercial  alloy  Rene  N4,  volume  fraction 
solid  as  a  function  of  height  docs  not  vary 
linearly  with  the  temperature  gradient  as  might 
be  expected.  This  strongly  affects  the 
permeability  as  a  function  of  depth  throughout 
the  mushy  zone  and  particularly  in  the  upper 
25%.  Using  the  3D  reconstructions  we  find  that  at  the  solid-liquid  solidification  front  the 
connectivity  of  the  uppermost  body  of  interdendritic  liquid  during  solidification  contains 
over  ninety  percent  of  the  total  interdendritic  liquid  and  extends  down  through  nearly  half  of 
the  entire  mushy  zone  height.  Using  these  reconstructions,  the  permeability  within  these 
structures  can  be  calculated  based  upon  the  observed  pressure  differential  resulting  from  fluid 
flow.  The  associated  permeability  values  calculated  offer  reasonable  agreement  with  similar 
treatments  of  permeability  in  the  literature.  We  find  that  cross-flow  transverse  to  the  primary 
dendrites  is  much  more  restricted  that  the  flow  parallel  to  the  growth  direction. 


Figure  3:  Streamlines  of  cross  flow 
through  a  volumetric  cross-section  of  the 
mushy  zone.  Solidification  direction  is 
normal  to  the  view.  Velocity  magnitudes 
are  presented  in  m/s. 


Composition  of  Mushes 

We  investigated  the  distribution  of  elements  within  freckles  and  in  their  surrounding  matrix 
in  Ni-based  supcralloys  applying  local-electrode  atom-probe  (LEAP)  tomography  combined 
with  scanning  and  transmission  electron  microscopy  techniques.  Two  directionally-solidified 
specimens  were  studied,  namely  ME-9  and  ME- 15  (S.  Tin  and  T.  M.  Pollock,  Metal.  Mater. 
Tram.  A  2003;  34A:1953).  These  specimens  contain  the  elements:  Ni,  Al,  Cr,  Co,  W,  Mo, 
Ta,  Rc,  Hf,  C,  and  B  have  the  nominal  compositions  specified  by  Tin  and  Pollock.  Both 
alloys  were  investigated  in  their  as-cast  state  to  preserve  the  original  partitioning  of  elements 
occurring  during  solidification. 

Both  misoriented  grains  and  freckles  were  present  in  directionally  solidified  ingots. 
Misoriented  grains  are  formed  in  the  liquid  close  to  the  solidification  front.  Therefore,  they 
arc  relatively  large  (on  the  same  order  of  magnitude  as  the  ingot  diameter),  uniformly- 
distributed.  and  have  a  composition  similar  to  the  average  value  in  the  solid.  The  mechanism 
of  formation  of  freckles  is,  however,  different.  Since  they  are  formed  by  thermosolutal 
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convection,  they  arc  smaller,  aligned  as  chains  close  to  an  ingot’s  perimeter,  and  their 
composition  is  close  to  that  of  the  last  liquid  remaining  prior  to  solidification:  the  latter  is  the 
same  as  of  the  intcr-dcndritic  regions  of  the  matrix. 

We  determined  the  identity  of  the  misoriented  defects  visually  as  well  as  by  compositional 
measurements  using  electron 
energy  dispersive  spectroscopy 
(EDS).  Using  EDS  we  compared 
the  composition  of  the  alloy  in  four 
different  regions:  dendritic  cores 
and  inter-dendritic  regions  of  both 
the  matrix  and  misoriented  defects 
(Fig.  4a).  We  utilized  the  four 
elements  W,  Re,  Ta,  and  Mo  as 
indicators  of  the  partitioning 
behavior  during  solidification 
because  it  is  known  that  Ta  and 
Mo  segregate  to  the  liquid  phase, 
whereas  W  and  Re  invcrsely- 
segregate  to  the  solid.  We  observed  both  types  of  partitioning  behavior.  In  the  freckles, 
similar  compositions  were  measured  in  regions  B,  C,  and  D,  which  arc  different  from  that  of 
region  A  (Fig.  4b).  The  average  composition  of  the  misoriented  grains  was,  how  ever,  about 
the  same  as  that  of  the  matrix  (Fig.  4c).  In  both  cases,  the  concentrations  of  W  and  Re  arc 
larger  in  the  dendritie  eores  than  in  the  inter-dendritic  region,  whereas  Ta  and  Mo  exhibit  the 
opposite  trend. 

The  microstrueture  of  both  alloys  is  not  homogeneous  at  all  pertinent  length  scales. 
Therefore,  we  decided  to  reconsider  the  previously-proposed  methodology  of  sample 
preparation  for  LEAP  tomography,  which  was  cutting  an  array  of  tips  from  a  freckle  and  its 
immediate  vicinity.  It  is,  however,  valuable  to  employ  the  conventional  method  (cutting  thin 
rods  and  electrochemically  sharpening)  for  precise 
measurements  of  the  average  composition  in  large 
misoriented  grains.  Both  the  matrix  and  a 
misoriented  grain  for  alloy  ME-9  were  analyzed; 

however,  results  were  eolleeted  only  for  the  yY  ’- 
phase  since  it  occupies  the  majority  of  the  volume. 

Only  one  tip  exhibited  a  mixture  of  both  the  y-Y 
and  yY  ’-phases,  from  which  we  draw  a  conclusion 
about  the  average  composition  of  the  grain.  The 
phases  can  be  determined  from  the  different 
partitioning  behavior  of  the  elements.  Fig.  5  is  an 
example  of  a  3-D  reconstruction,  showing  the 
distribution  of  three  atomic  species,  for  the  sake  of 
elarity:  A1  (red  dots),  Cr  (blue  dots),  and  C 

(spheres).  The  interface  between  the  y-Y  and  y Y  ’- 
phase  was  determined  using  a  9.3  at.%  Cr 


Figure  5:  A  3-D  atom-probe 
tomographic  (APT)  reconstruction 
showing  the  y-  and  y’-phases.  The 
phases  are  identified  from  the  high 
partitioning  of  Cr  (blue)  to  the  y-,  and 
of  A1  (red)  to  the  y’-phase.  C 
partitions  to  the  y-phase.  Both  phases 
are  separated  by  a  9.3  at.%  Cr  iso- 


6 


isoconcentration  surface.  The  other  elements  (not  presented  herein)  exhibit  different 
partitioning  behaviors:  Cr,  Co,  W,  Re,  Mo,  Hf,  B,  and  C  partition  to  the  y-phase,  while  Ni, 
Al,  and  Ta  -  to  the  y’-phasc.  The  partitioning  ratios  of  each  of  the  elements  were  also 
determined.  Although  this  trend  is  general,  there  are  some  exceptions  in  which  the 
partitioning  behavior  of  minority  elements  is  reversed.  For  example,  we  find  the  (slight) 
preference  of  W  for  the  y-phase  is  reversed  in  favor  of  the  y’-phase  as  a  result  of  the  presence 
of  Ta.  We  hypothesized  that  because  of  tantalum’s  preference  for  the  Al-sublattice  sites  of 
the  y’-phase,  Ta  has  a  stronger  tendency  to  occupy  these  sites,  thus  rejecting  the  W  atoms  to 
the  y-phase.  We  confirmed  this  hypothesis  employing  first-principles  calculations,  using 
VASP,  of  the  substitutional  formation  energies  of  W  and  Ta  at  different  sites  in  both  the  y- 
and  y’-phases.  The  calculations  demonstrated  that  the  energetic  driving  force  for  Ta 
partitioning  to  the  y’-  phase  is  twice  that  of  W. 


The  additional  method  we  utilized  to  investigate  the  composition  in  the  four  regions  of 
interest  (see  Fig.  4)  involves  lifting-out  samples  from  these  regions  using  a  FEI  Helios  dual¬ 
beam  focused-ion  beam  (FIB)  microscope  (FIB).  This  procedure  is  illustrated  in  Fig. 6. 


In  this  framework, 
we  developed  a 
method  that 

enables  us  to  take 
advantage  of  the 
superior 

detectability  and 
mass  resolution  of 
the  LEAP 

tomograph,  and 
the  ability  to 
quantify  elemental 
concentrations  at 
much  longer 

length  scales  (10- 
100  pm)  than 
those  typically 
obtained  using 
atom-probe 
tomography  ( 1 0- 
100  nm). 


Figure  6:  Producing  atom-probe  tomographic  microtips  on  a  Si  micropost 
array,  (a)  A  dendrite  core  in  a  freckle,  from  which  a  wedge  is  taken,  is 
marked  by  a  circle.  The  subsequent  step  involves  cutting  a  wedge  and 
attaching  it  to  a  micro-manipulator  that  lifts  it  out  of  the  bulk,  (b)  The 
wedge  is  attached  to  a  pre-sharpened  silicon  microtip,  and  then  sharpened 
by  ion-milling  to  obtain  a  40  nm  radius  of  curvature,  which  is  required  for 
LEAP  tomography. 


Our  technique  is  based  on  a  mass  balance  existing  between  the  y-  and  y’-phases  (Eq.  1 ),  and  is 
applicable  for  all  LEAP  3-D  reconstructions  comprising  both  the  y-  and  y’-phascs,  Fig.  5. 
The  difference  between  the  average  concentration  of  an  element  /  in  the  dendrite  cores  (Cavg  - 
measured  by  EDS)  and  its  concentration  in  the  y-phase  alone  ( Cy ,  measured  by  APT)  is 
plotted  against  the  difference  between  the  content  of  i  in  y ’(Cy  )  and  y,  respectively. 


C„-Cr+(p,-Cr>*,  (1) 
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This  procedure  is  made  for  all  elements  detected  employing  EDS.  The  result  is  a  linear 
dependence  with  a  slope  equals  to  the  volume  fraction  of  y '(<p  ).  A  very  useful  application  of 

this  seheme,  besides  indieating  the  consistency  between  EDS  and  APT,  is  that,  in  this 
manner,  we  ean  derive  the  average  composition  of  any  element  in  the  dendrite  eore  (or  any 
of  the  other  4  regions  of  interest),  even  for  those  not  detected  by  EDS:  see  results  in  Fig.  7. 

Based  on  the  above  methodology,  we  estimated  the  solid/liquid  partitioning  behavior  of  Hf. 
In  contrast  with  other  elements  (such  as  W,  Rc,  Ta,  Mo),  of  which  the  partitioning  behavior 
(and,  therefore,  their  role  in  freckle  formation)  is  known,  there  are  insufficient  data  regarding 
the  solid/liquid  partitioning  behavior  of  Hf,  since  the  Hf  content  is  usually  too  small  to  be 
quantified  by  conventional  means,  such 
as  EDS.  Although  the  nominal  Hf 
concentration  is  relatively  low  (0.05 
at.%),  the  question  of  Hf  miero- 
segregation  to  the  solid/liquid  interface 
during  solidification  is  important,  since  it 
is  a  heavy  element  (178.49  g  mole'1)  and 
is,  therefore,  more  likely  to  eause  density 
inversion  that  results  in  thermosolutal 
eonveetion  and  freckle  formation.  The 
result  is  an  average  concentration  of  Hf 
in  the  matrix  dendrite  eores  of  0.10±0.01 
at.  %  and  0.148±0.012  at.  %  for  the  ME- 
9  and  ME- 15  alloys,  respectively.  These 
values  are  mueh  larger  than  the  nominal 
Hf  concentration  of  both  samples,  0.05 
at.%,  indieating  that  Hf  partitions  to  the 
dendrite  eores  (solid).  Moreover,  no  Hf 
was  detected  in  the  samples  lifted-out 
from  the  freckles  dendrite  cores.  We 
eonelude  that  Hf,  by  partitioning  to  the  solid  during  solidification,  reduees  the  average 
density  in  the  liquid  and  this  may  enhance  freekle  formation.  Quantitatively,  we  evaluated 
the  liquid  density  using  a  semi-empirieal  expression  for  the  density  of  a  multi-eomponent 
liquid  (P.  K.  Sung,  D.  R.  Poirier,  and  E.  McBride,  Mater.  Sci.  Eng.  A  231  (1-2):  189,  1997), 
whieh  takes  into  aeeount  a  mixing  term,  A  V*4,  Eq.(2): 

The  liquid  density  decreases  by  0.12%  due  to  Hf  depletion,  whereas  observations  (R.A. 
Hobbs,  S.  Tin,  C.M.F.  Rac,  Metall.  Mater.  Trans.  36A:  2761,  2005)  indicate  that  freckles 
appear  only  when  the  density  inversion  term  exeeeds  2-5%. 

Using  the  above  techniques  we  quantified  the  compositions  in  the  dendrite  eores  of  the 
matrix  and  freekles  of  the  ME- 15  alloy  using  LEAP  tomography.  The  other  regions  (inter- 
dendritie),  were  analyzed  only  using  EDS.  Based  on  Eq.(2),  we  estimated  the  density 
inversion  term,  (Ap/p0),  to  be  2%  and  0.5%  in  the  Me- 15  and  ME-9  alloys,  respectively. 
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Figure  7:  A  mass-balance  diagram  plotted  for  the 
matrix  dendrite  cores  in  the  ME-9  alloy.  The 
elements  shown  are  detected  by  EDS,  while  the 
concentration  of  the  other  minority  elements  can  be 
interpolated  using  the  fitting  formula. 
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Theory  of  Flow  in  Mushy  Zones 


We  examined  the  convective  state  due  to  weak  sidewall  heat  losses  in  a  mushy  zone  during 
the  steady  directional  solidification  of  a  binary  alloy.  The  configuration  consisted  of  a  warm 
liquid  region  and  a  cold  solid  region  separated  by  a  mixed  phase  region,  the  mushy  zone, 
which  is  modeled  as  a  reactive  porous  matrix.  The  structure  of  the  convection  that  arises 
from  horizontal  temperature  gradients,  induced  by  the  heat  losses  at  the  sidewalls,  is 
characterized  by  a  set  of  nondimensional  parameters  that  describe  the  effects  of  latent  heat, 
composition,  permeability,  and  thermal  and  solutal  buoyancy.  We  observed  a  wide  range  of 
behaviors  and  showed  that,  as  the  critical  Rayleigh  number  for  convection  in  a  horizontally 
uniform  mush  is  reached,  we  begin  to  see  the  precursors  of  chimneys  near  the  cooled 
boundaries.  The  strength  of  the  cooling  plays  an  important  role  in  determining  the  strength 
and  degree  of  localization  of  the  convection  near  the  boundary.  We  find,  in  common  with 
other  authors,  that  upflow,  in  this  case  caused  by  lighter  fluid  being  released  at  the  cooled 
sidewalls,  leads  to  regions  of  dissolution,  which  are  precursors  to  chimney  formation.  We 
also  indentified  the  effects  of  the  different  physical  parameters  on  the  steady  states 
convection  patterns. 

The  above  mentioned  stability  theory  shows  that,  indeed,  nonlinear  instabilities  grow 
preferentially  near  the  walls,  where  the  up  flow  melts  away  the  dendritic  structure  producing 
channels.  The  conditions  for  chimney  growth  depend  not  only  on  the  local  Rayleigh  Number 
the  measures  buoyancy,  but  the  Biot  Number  that  measure  the  heat  losses. 

Most  importantly,  this  paper  showed  why  a  correlation  depending  on  Rayleigh  Number  only 
reasonably  predicts  the  conditions  of  chimney  growth  when  the  front  is  flat,  but  fails  when  it 
is  curved  (here,  by  heat  losses).  The  proper  criterion  depends  both  on  Rayleigh  Number  and  a 
second  parameter  that  measures  front  curvature  and  hence  horizontal  convection. 

We  also  examined  the  dynamics  of  a  mushy  layer  in  directional  solidification  for  the  case  of 
a  thin,  near-eutectic  mush  with  a  deformable  and  permeable  mush-liquid  interface.  We 
focused  on  the  onset  of  convection  using  linear  stability  analysis,  and  the  weakly  nonlinear 
growth  of  liquid  inclusions  that  signal  the  onset  of  chimneys.  This  analysis  was  compared  to 
past  analyses  in  which  the  mush-liquid  interface  is  replaced  by  a  rigid,  impermeable  lid.  We 
found  qualitative  agreement  between  the  two  models  but  the  rigid-lid  approximation  gave 
substantially  different  quantitative  behavior.  In  linear  theory  the  rigid-lid  approximation 
leads  to  an  over-estimate  of  the  critical  Rayleigh  number  and  wave  number  of  the  instability 
by  as  much  as  a  factor  of  about  5  and  2  respectively.  The  onset  of  oscillatory  instability  was 
also  changed  by  a  factor  of  about  5  in  composition  number  C.  In  the  weakly-nonlinear 
theory,  the  location  of  the  onset  of  liquid  inclusions  is  near  the  undisturbed  front  for  the  free¬ 
boundary  analysis  whereas  it  lies  at  the  centre  of  the  mushy  layer  when  the  rigid-lid 
approximation  is  used.  An  example  of  the  onset  of  chimneys  that  follow  from  this  calculation 
is  given  in  Figure  8. 
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